Full-Length Genome Sequence of a Novel European Antigenic Variant Strain of Infectious Bursal Disease Virus

ABSTRACT We report the full-length genome sequence (compared to reference sequences) of a novel European variant strain of infectious bursal disease virus (IBDV), designated 19P009381 (AxB1). This should help to further identify such viruses in Europe.

I nfectious bursal disease (IBD) is a ubiquitous immunosuppressive and sometimes lethal disease of young chickens that is induced by IBD virus (IBDV), a member of the Birnaviridae family and Avibirnavirus genus. IBDV has two double-stranded segments of RNA (1) (A and B). Segment A encodes, among others, the capsid protein VP2. VP2 contains a hypervariable region (HVR), itself containing four loops (2). Segment B encodes the viral polymerase VP1. Viruses from genogroup A2B1 (3,4) are known as antigenic variants. Sequence information on these strains in Europe is scarce (5,6).
In 2019 in Brittany, France, five bursae of Fabricius (BFs) with lesions evocative of IBDV were pooled for viral RNA extraction using the QIAmp viral RNA kit (Qiagen) and deep sequenced using Ion Proton technology, as described previously, with the Ion Total transcriptome sequencing (high-throughput RNA sequencing) kit (7). A total of 2,387,143 reads, with an average length of 95 nucleotides, were obtained. All tools were run with default parameters unless otherwise specified. Quality assessment was performed with fastqc v0.11.9 and trimming with Trimmomatic v0.36 (options: ILLUMINACLIP=oligo_file.fasta:2:30:5:1:true, LEADING=3, TRAILING=3, MAXINFO=40:0.2, and MINLEN=36). Reference-based assembly using BWA-MEM v0.7.8 and SPAdes v3.10.0 for a de novo assembly step and then alignment of the de novo consensus sequences to references (GenBank accession number X84034.1 for segment A and GenBank accession number AY459321.1 for segment B) yielded median coverage values of 232Â and 154Â for segment A (3,261 bp, with a GC content of 53.3%) and segment B (2,827 bp, with a GC content of 53%), respectively. Coverage in the 59 and 39 untranslated regions (UTRs) was lower than that in the rest of the genome for both segments (7Â and 45Â for segment A extremities and 11Â and 14Â for segment B extremities, respectively), and thus the sequences were obtained by alignment to references. The full-genome consensus sequence (reference number 19P009381) was corrected based on actual reads, using the Integrated Genomic Viewer (IGV) v2. 8 1, Algeria, 2015), with 97.40% and 97.18% nucleotide identity, respectively. Phylogenetic analyses showed that segment A does not cluster robustly with any described genogroup (Fig. 1), including the A2-US antigenic variant (evolutionary distance of 0.088), although it appears more closely related to a small set of segment A sequences from The evolutionary history was inferred using the neighbor-joining method (9). The percentages of replicate trees in which the associated taxa clustered together in the bootstrap test (1,000 replicates) are shown next to the branches (10). The evolutionary distances were computed using the Kimura twoparameter method (11) and are in the units of the number of base substitutions per site. The analysis involved 492-and 461-nucleotide sequences for segment A and segment B, respectively. All positions containing gaps and missing data were eliminated. There were totals of 355 and 461 positions in the final data sets for segments A and B, respectively. Evolutionary analyses were conducted in MEGA7 (12).
Portugal. However, segment A of strain 19P009381 diverges from previously defined A lineages by less than 9%. Thus, although it appears to represent genetic elements indicating a new genogroup, more sequences are required for confirmation. Segment B of this virus clustered with genogroup B1 (Fig. 1).
VP2 HVR alignments of strain 19P009381 with classic strains and documented antigenic variants confirmed the presence of mutations associated with altered antigenicity (Table 1).
In summary, BLASTn results, amino acid alignments, and phylogenetics support the assertion that French strain 19P009381 is a new example of a European antigenic variant, tentatively assigned to the AxB1 genogroup.
Data availability. Raw next-generation sequencing (NGS) reads and full-length and partial sequences for segments A and B were deposited in the SRA and GenBank under accession numbers SRR15808265, OK043826, OK043827, MZ556959, and MZ556960.